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1^ ■ 1 INTRODUCTION 



Conventional theo ry of accretion discs proposed by Shakura 



Sunyaev (1973) was successful in interpretation of obser- 
vations of the accretion discs in the binary systems. This 
theory satisfactory predicts the general properties of the 
discs around compact objects. Nevertheless some important 
phenomena connected with the accretion appeared incom- 
patible with the conventional theory. The most important 
among them are the accretion discs with anomalously low 
luminosity and jet-like outflows from the objects. The best 
example in this regard is AGN M87. 

Advection dominated disc model (ADAF model ) was 
proposed to explain the under luminous discs, i n particular 
dis c around BH located in t he Galactic Centre (Narayan & 





1995 



solve the problem of jet-like outflow especially in the cases 
when the kinetic luminosity of jets is comparable (like in 
the case of SS433 case) or even exceeds ( like in the case of 
the jet from M87) the bolometric luminosity of the object. 
In fact, these problems including low luminosity accretion 
discs, very efficient outflow from the discs and its collimation 
into jets could be internally connected. 

Observations show that jets are strictly connected with 
the disc accretion. In all jet detections a signature of a disc 



ABSTRACT 

We consider an extreme case of disc accretion onto a gravitating centre when the 

viscosity in the disc is neghgible. The angular momentum and the rotational energy of 

the accreted matter is carried out by a magnetized wind outflowing from the disc. The 

outflow of matter from the disc occurs due to the Blandford & Payne(1982) centrifugal 

mechanism. The disc is assumed to be cold. Accretion and outflow are connected by 

the conservation of the energy, mass and the angular momentum. The basic properties 

of the outflow, angular momentum flux and energy flux per particle in the wind, do not 

depend on the details of the structure of the accretion disc. In the case of selfsimilar 

accretion/outflow, the dependence of the rate of accretion M in the disc depends on 

1 

the disc radius r on the law M ~ r ^("^-i) , where a is a dimensionless Alfvenic radius. 
In the case of a 3> 1, the accretion in the disc is provided by very weak matter outflow 
from the disc and the outflow predominantly occurs from the very central part of the 
disc. The solution obtained in the work provides mechanism which transforms the 
gravitational energy of the accreted matter into the energy of the outflowing wind 
with efflciency close to 100%. The final velocity can essentially exceed Kepler velocity 
at the site of the wind launch. This mechanism allows us to understand the nature of 
the astrophysical objects with low luminosity discs and energetic jet-like outflows. 
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accretion was found as well. It is important for understand- 
ing of the mechanism of the jet ejection that this phenomena 
is not connected with the specific nature of the central ob- 
ject. Jets are formed irrespectively to the fact that the cen- 
tral object is black hole ( like in the case of jet from AGNs), 
neutron stars or protostar as it takes place in all cases of 
the jets from Young stellar objects. It is reasonable to as- 
sume that the jets are directly connected with the accretion 
mechanism itself rather than with the nature of the central 
object. However, in this case one observational fact needs to 
be explained. In all cases when it was possible to observe the 
base of the jets it appears that the jets are launched from 
the very central part of the disc which produce impression 
that central object could be connected some way with the 
process of the jet ejection. At least the fact that jet is ejected 
from the very central part of the accretion disc rather then 
from all the disc surface demands explanation. 



Narayan et al. 1998). This theory, however does not 



Another difficult problem is that in some cases the 
mechanism of ejection is surprisingly efficient. For exam- 
ple, total bolometric luminosity of M87 does not exceed 



10 ergs/s (Biretta Sern & Harris 1991), while the total ki- 
netic luminosity of the jet from M87 is as high as lO^'^ergs/s 



(Bicknell & Begelman 



199e 



Reynolds et al. 1996). Thus, if 



to estimate the gravita tional energy release on the basis of 
conventional theory of ghakura fc Sunyaev ( [l973|) , the ki- 
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netic energy luminosity of the jet from M87 exceeds by two 
orders of magnitudes the gravitational energy released at 
the accretion. The conventional models of the disc accretion 
(Shakura & Sunyaev type or ADAFs) do not predict the 
existence of such objects. The example with M87 jet shows 
that one needs to explore new regimes of accretion. 

The most evident modification is incorporation into the 
model of the magnetic field. It has already been widely rec- 
ognized that the magnetic field has important impact and 
may play leading role in the ejection of the plasma from 
accretion discs. In this regard, two pro cesses are of special 
interest. First one is related the idea of Blandford & Payne 



(1982). They have demonstrated that the magnetic field re- 
sults into instability of the particles at the Kepler orbit. If 
the angle between the force line of the magnetic field and 
the disc plane is less than 60° the particles are freely ejected 
from the dis c by centrifugal for c e. Th e second process was 
proposed by Pelletier & Pudritz (1992). They argue that the 



winds from the accretion disc can carry out noticeable part 
of the angular momentum of the accreted material increas- 
ing the accretion rate. Together, these work clearly demon- 
strate that the magnetic field of the disc results into the 
outflow from the disc at rather general conditions and this 
wind can carry out essential part of the accreted material 
angular momentum. 

In this work we consider an extreme case when the an- 
gular momentum of the disc is carried out by the wind. The 
viscous stresses are fully neglected. We solve the problem 
of the disc accretion under these conditions selfconsistently 
with the problem of the wind outflow from the disc. Fortu- 
nately, to provide selfconsistency of the processes of accre- 
tion and outflow one do not need to know information about 
the detailed structure of the accretion disc. The laws of con- 
servation of mass, energy and angular momentum appeared 
sufficient to provide the selfconsistency of the processes of 
accretion and outflow. 

The paper is organized as follows. In the flrst section 
we discuss the qual i tative picture of the outflow due to the 
Blandford fc Payne (1982) centrifugal instability and related 
structure of the poloidal magnetic field. In the second sec- 
tion the basic equations and connection of the accretion and 
outflow are discussed based on of the conservation laws. In 
the third section the selfsimilar solutions of the problem are 
presented. In the last section we discuss the physical sense 
of the solution. 



2 QUALITATIVE PICTURE OF THE 

ACCRETION UNDER AFFECT OF WIND. 

The idea that the wind can carry out essential fraction of 
t he angular momentum fr o m the disc has been explo red in 



( [Ferreira fc Pelletiei| pjg^ ; [Combet fc Ferreira |pOO^ ). The 



approach of Ferreira is analog o f the pioneering studies by 
Bisnovatyi-Kogan & Blinnikov (1977) of the disc accretion 



onto gravitating centre from a magnetized cloud surrounding 
the gravitating centre. They considered the accretion as the 
process of diffusion of the matter across the magnetic fleld 
lines rooted into interstellar medium. Actually, t his kind of 
flow is similar to the Hartman flow in a channel Landau & 




Figure 1. The structure of the magnetic field in the accretion disc 
and in the outflowing wind. The disc is shadowed. The magnetic 
field lines in the disc are distributed chaotically. At the base of 
the wind from the disc all the magnetic field lines are opened. 
Their direction is accidental as well, therefore the total magnetic 
flux leaving one side of the disc equals to zero. The box drown in 
black thick lines is the region of integrations of conservation laws 
connecting the properties of the disc and the wind. 



assume in this scenario an existence of a mechanism which 
reduces the electric conductivity of plasma on a few orders 
of magnitude to provide efficient accretion. This scenario 
has another problem. The velocity of the plasma outflow 
from the disc exceeds the velocity of all MHD perturba- 
tions. Therefore no signal from the cloud can propagate to 
the disc along the fleld lines. Therefore, there is an internal 
inconsistency of the model. From one side it requires that 
the fleld lines are rooted into the cloud to prevent advec- 
tion of the fleld lines with the matter. On the other hand 
the cloud can not affect on the magnetic fleld in the disc 
because it is causally disconnected from it. 

Here we suggest another approach which does not face 
these kind of problems. The source of the magnetic field in 
our model is the plasma of the accretion disc itself like Sun 
is the source of the magnetic field on the solar surface and 
in the interplanetary space. It is reasonable to expect that 
in the disc the magnetic fi eld lines are distrib u ted r ather 
chaotically as suggested in Blandford & Payne (1982) and 
is shown in Fig. hi. 

In the physical model of Ferreira et al. the matter have 
to diffuse across the poloidal magnetic field lines. In our 
approach the matter falls down onto the gravitation centre 
together with the frozen-in magnetic field. Since dissipation 
processes are not necessary for this, the process of accretion 
and outfiow can be considered in the limit of ideal magneto- 
hydrodynamics where no heating of the matter. Therefore, 
the disc is cooled to low temperature even if the plasma was 
hot initially. 

Thus, it is reasonably to expect that the disc is geomet- 
rically thin with thickness 



(1) 



Lifshitz (1982). Taking into account that the plasma typi- 
cally has very high electric conductivity it is necessary to 



where Ca is the sound velocity in the disc, Vk is the Kepler 
velocity of the disc and r is the cylindrical radius. The out- 
fiow of plasma from the cold disc occurs due to mechanism 
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specified by Blandford & Payne (1982). If the angle between 
the magnetic field line and the plane of the disc is less than 
60° then a particle motion on the Kepler orbit becomes un- 
stable. Any perturbation of the radius of the orbit results 
into sliding of the particle along the field line outward from 
the gravitating centre. Therefore, the outflow of the plasma 
from the accretion disc can occur even if the disc is cold. 



3 EQUATIONS. 
3.1 Basic equations 

In this work we study a steady state axisymmetric flow of 
an ideal plasma. The condition of ideality has a form 



E+ -V X B = 0, 



(2) 



where E is the electric field, B - magnetic field and v is the 
velocity of the plasma. 

It is convenient to consider MHD equations for ideal 
plasma in the form o f conservation laws. According to 



Lan- 



dau & Lifshitz (L988) the energy- momentum of plasma sat- 



isfies to the equation 



0, 



where Tik is the energy momentum tensor in the form 



(3) 



(4) 



r,2 

Where Too = e + is the sum of the thermal and magnetic 
field energy densities. We consider here only nonrelativistic 
flows. Therefore the term in Too connected with the electric 
field is omitted. Spatial components of the tensor aik are as 
follows 



Too 




92 


93 


gi 




0"12 


0-13 




0"21 


0"22 


0"23 


93 


0"31 


0"32 


0"33 



CTife = pViVk + pSik 



— (BiBk 
47r^ 



(5) 



Equation (^) still does not contain gravitational force. 
To take into account the gravitational force it is sufficient to 
add into right hand part of Eq. Q with a spatial components 
-pGM-S- so that the equation for atk 



of the tensor a term 
becomes as follows 



daik 



^pGM 



(6) 



Here Ri is the component of radius- vector from the position 
of the gravitating centre. 

The energy density fiux qi is as follows 



(7) 



In this expression the gravitational field is taken into account 
explicitly. Therefore the equation for the energy density flux 
remains as follows 

1^=0 (8) 
axk 

for the steady state flow. 

In the Cartesian coordinates the angu lar momentum 
tensor is i ntroduced as rriik — eimpXmTpk (Landau & Lif- 
shitz 1988 ), where Simp is unit antisymmetric tensor. In this 



case the density flux of the z component of the angular mo- 
mentum is as follows 



k — pvivv^p — —rBiBu: 



(9) 



The conservation equations for the matter and the mag- 
netic fluxes are 
dpvk 



= 0, 



(10) 



= 0. 



(11) 



dxk 
and 
dBk 
dxk 

These equations are supplemented by the condition of 
stationarity of the magnetic field 

curlE = 0. (12) 
This condition means that ^ — 0. 



3.2 The role of the induction equation in the 
problem of accretion. 

In the limit of axisymmetric flow an azimuthal component of 
Eq.dlil) and frozen- in condition (^ for the same component 
give a couple of equations 

d_ 

dr 
and 

E^ + - (v^Br - VrB^) = 0. 

c 



{rE^) = 0, 



(13) 



(14) 



Here r is the distance from the symmetry axis. 
The solution of Eq.([3|) gives that 



(15) 



where A is some constant. This solution diverges at r 0. 
On this reason it is assumed that Eip = (Mestel 1999 ). This 
is a conventional assumption for study of axisymmetrical 
MHD winds. 



Contopoulos ( 1996 ) was the first who pointed out that 
if we deal with the accretion of an ideal plasma onto a grav- 
itating centre E^ is not equal to zero. Indeed, as it follows 
from Eq.([l4|), E^ at the equator is equal to ^VrBz because 
Br ~ and = Q a.t the centre of the disc. At accretion 
Vr 7^ 0, therefore, E^ 7^ as well. Thus, in the region of the 
accretion fiow Etp can not be neglected because it is con- 
nected directly with the radial velocity of the plasma in the 
disc. 

At first glance the advection of the magnetic field to the 
gravitating centre results into accumulation of the magnetic 
field lines at the centre and to infinite growth of the magnetic 
field. This happens if the direction of the magnetic field is the 
same at one side of the disc as it is shown in Fig. |2| Actually 
this is unrealistic case. The direction of the magnetic field 
on one side of the disc varies accidentally so that in average 
the magnetic flux leaving one side of the disc is zero. That is 
why at the advection of this magnetic fleld to the gravitating 
centre there is no accumulation of the magnetic field fiux 
because it equals to zero in average on the surface of the 
disc. The infinite growth of the magnetic field at the centre 
is also prevented by the reconnection of the magnetic field 
lines as it shown in Fig. ^ 



4 S.V. Bogovalov and S.R. Kelner 




to take into account that w ~ 14 and that vi 



the 



Figure 2. The advection of the magnetic field lines of one direc- 
tion to the rotational axis results into infinite accumulation of the 
magnetic field flux. There is no way to annihilate the magnetic 
field with this topology. 




Figure 3. The advection of the magnetic field lines of differ- 
ent directions to the rotational axis does not result into infinite 
accumulation of the magnetic field flux because in average the 
magnetic flux on one side of the disc equals to zero and there is 
way to annihilate the field lines in the process of reconnection. 



As was noted in Bogovalov fc Tsingano!; ( 2000 ) , the dy 



namics of an ideal plasma is invariant with respect to the 
redirection of the magnetic field lines. Therefore the dynam- 
ics of plasma in the fields shown in Fig. |^ and Fig. ^is the 
same. With this in mind, we are able to consider an out- 
fiow of plasma from the accretion disc with the azimuthally 
symmetrical magnetic field shown schematically in Fig. e1 



3.3 Along field line MHD equations of the wind 

We consider the case when E^p ^ 0. As it follows from 
Eq.(p^ in this case the velocity has a poloidal component 
v± orthogonal to the poloidal magnetic field line. Neverthe- 
less, this component of the electric field can be neglected at 
the consideration of dynamics of the wind under the condi- 
tion v_i_ <^ V, where v is the full velocity of the plasma. If 



toroidal electric field can be neglected under the condition 



BpT- 



14 



< 1. 



(16) 



It foUows from this condition that the toroidal elec- 
tric field can be neglected if the radial velocity in the disc 
Vr <^ Vk- But this is not enough to provide the condition 
(p^. The product BpV should not drop down too strongly 
compared with the initial value. This means that if for exam- 
ple, poloidal field falls down as there is limitation on the 
distance where the toroidal electric field can be neglected. 
Further we will be interested in the solution of the problem 
at the distance small compared with the size of the accretion 
disc. In this region can be neglected. The conclusion that 
in the limit Vr <C Vk affect of the toroidal magnetic field on 
dynamics of the wind can be neglected has been shown ex- 
plicitl y for the selfsimilar solutions by Bogovalov & Kelner 

( poosj ). 

In the region where E^p — the poloidal velocity is di- 
rected along the poloidal magnetic field. This directly follows 
from Eq.(|l4|). In this case we have 



Ip = prVp Vp 



Br, 



Anpvp 



(17) 



for the angular momentum flux density along a poloidal flled 
line and 



1p = P^'P I Y 



GM 



- Qr 



Br, 



Anpvp 



(18) 



for the energy density flux along a poloidal fleld line. If to 
take into account that the fluxes of the angular momentum 
(JtppdS, energy QpdS, matter rpvp and magnetic field fiux 
BpdS are conserved as it is demonstrated in Fig. ^, it can 
be obtained that the following two integrals of motion take 
place along the field lines 



where / = AupVp/ Bp and 
11^ GM 1 



R 



■QrB^ 



E. 



(19) 



(20) 



The first equation from this couple is the conservation of the 
angular momentum per particle and the second one is the 
conservation of the energy per particle along a field line. 

The frozen-in condition for the poloidal component of 
the electric field gives that 

rQBp + VpBip = VipBp. (21) 

We take into account in Eq.(^) that due to conservation 
of the magnetic field fiux BpdS and product Epdl, which is 
constant due to induction equation (h2) there are relation- 
ships Er = —rQB^/c and Ez = rD.Br/c. Combining Eq.dil]) 
with Eq.(|l9|) gives that 



Lfvp — r^QBp 
fvp - Bp 



(22) 



Denominator of this expression goes to zero at the Alfvenic 
point where Vp = Bp/^/^np. The nominator of this expres- 
sion must equal to in this point to provide regularity of 
Vip. From this condition we obtain that the momentum per 
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Figure 4. Fluxes of the magnetic field, matter, energy and an- 
gular momentum between any two close field lines with the cross 
section dS = 2-nrdl are conserved. The condition curl E = gives 
that the product Epdl, where dl is the distance between these 
field lines is also conserved. 



particle equals to L 
Alfvenic point. 



D, r. , where r a is the radius at the 



3.4 The connection of the accretion and outflow 

Irrespective to the inner structure of the disc, the disc and 
the wind are connected by the conservation laws. Let us 
consider a fragment of the disc as it is shown in Fig. |l| 
A rectangular region is shown by thick black lines. The upper 
and lower boundaries of this region are located at the base 
of the wind from the disc. All the magnetic field lines of the 
wind are rooted here. 

Conservation laws in integral form are 



Ik dSk = 0, 
for the angular momentum, and 
pvk dSk — 



(23) 



(24) 



for the mass conservation. Integration here is performed 
along a closed surface surrounding the volume. 

Let us perform integration of the angular momentum 
flux over a surface which includes a small part of the disc 
shown in Fig. ^. The surface consists of the sides Sl at radius 
ri, side Sr at radius r2 and upper and down sides Sup and 
Sd- Integration gives 



h/2 



r pVtpVr dz 



h/2 



+ 



+ 



47r 



h/2 



h/2 



h/2 
-h/2 

4-K 



r pVtpVr dz 



+ 



r BrBtp dz 



+ 



(25) 



+2 (^rpv^Vz — -^rB^Bz 



rdr — 0. 



Integration across the disc is performed in the interval on 



z from 



to 



correspondingly at the radiuses rl and 



r2. Integrations along Sup and Sd are equal to each other 
because the vector dS and the component of the velocity 
Vz change sign simultaneously. Therefore, we simply double 
the integration along the surface Sup- The terms containing 

1 rh/2 o 

magnetic field of the disc ^ j^h/2 ^ ^rB^ dz are much less 

r-h/2 2 

compared with the terms J_^y2 r pv^Vr dz provided that the 
condition 

pVi » -^B^ (26) 

takes place. In other words the angular momentum flux is 
fully dominated by matter rather than the magnetic field 
in the disc. Under this condition Eq.(p5|) is reduced to the 
differential form as follows 



d 
rdr 



+2 ( rpv^Vz 



h/2 



pv, 



dz 



h/2 



(27) 



— rB^B, 
An 



= 0. 



The subscripts disc and winds denote the variables describ- 
ing the disc and the wind at the base ( at the surface Sup)- 
According to this equation the angular momentum of the 
disc is carried out by the outflowing plasma and by the mag- 
netic stresses in the outfiow. 

After similar manipulations it is easy to obtain that the 
equation for r-component of the momentum is as follows 



r dr\ 



+ 2 pVrVz 



h/2 



pv^az J 



h/2 



- — BrB^ 
4% 



h/2 



pv^ dz 



h/2 



h/2 



pGM 



dz 



h/2 



(28) 



Here we neglect the thermal pressure, assuming that the disc 
is cold. This equation describes the balance of forces along 
radius in the disc. In the limit — » and under condition 
(bq) the Keplerian rotation takes place 



h/2 



pv^ 



, dz = 



1 



h/2 



h/2 



pdzGM- 



(29) 



h/2 



Here we also neglected the term pVrV^ in the wind. Below 
we will see that this term is also equal to zero because in the 
case under consideration poloidal velocity of the wind equals 
to zero at the base. Everywhere below we accept that the 
azimuthal velocity of the plasma in the disc Vtp = Vk- The 
conservation law for the energy flux is as follows 



qk dSk = 



(30) 



Integration of the energy flux over the surface shown in 
Fig. hi gives the following equation 



(31) 



+2 pv 



■ 0. 



The first term in this equation is the variation of the en- 
ergy flux in the disc. Under condition (^ the energy flux 
of the disc consists only on the kinetic and gravitational en- 
ergy of the matter. The magnetic field energy fiux can be 
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neglected. The energy from the disc is carried out by the 
wind. Radiation is neglected because the disc is cold. 

The last equation necessary to connect the accretion 
with outflow is the matter conservation 



d 
rdr ' 



[t pvrdz] +(2p^;,)_^ = 0. 

\ J_h/2 /disc 



(32) 



ft is convenient to introduce the accretion rate in the 

h/2 



disc 

M = — 27rr / pVrdz. 

J -h/2 

In this case the last Eq.(^2|) takes the form 
^■wrpVz = 0. 



dM 
dr 



If to insert this equation into (p^), we obtain 



or ^ ' 



dM 
dr 



r_B^ 
f 



= 0. 



(33) 



(34) 



(35) 



But the expression [rVk 1 equals to L = rj^Qk - 

V / wind 

the angular momentum per particle in the wind. Differential 
equation for M takes the form 



d ■ dM 9 

^irV,M)-^rArrn,{r)^0, 



(36) 



where (r) is the Alfven radius of the force line rooted into 
disc at the point with radius r, 



GM 

3~ 



(37) 



is the Keplerian angular velocity of the disc at the radius r. 
The solution of the equation (pm gives 



In 



M 



dr^ 



4(r^(r) — r^) ' 



(38) 



where Mmax is the accretion rate at the largest radius of the 
disc Tinax- Lct US cousider the simplest case when rA{r) = 
ar, where a is some constant. The physical meaning of a is 
simple. This is the Alfvenic radius of a force line expressed 
in the radius of the point at the disc where the force line is 
rooted. In this case the accretion rate in the disc varies with 
r as follows 



M = Mrr 



\ I max / 



(39) 



The accretion rate in the case under consideration varies 
with r according to power law. It is interesting that in the 
case a S> 1 the power index appears close to zero. This 
results into very interesting astrophysical implications which 
will be discussed below. Now we consider the role of Eq. (^l[) . 

Using (^^ it is easy to obtain from (^l|) that 



d_ MVj 
dr- 2 



d: 



GM VkB^ 



f 



= 0. (40) 



1 [-BXS]; 



We use notations that . — 

47r pv^ 

previous section and the condition 

CM ^hrBips 



introduced in the 
The term (^ - 



E. We are interested in the solutions 



which allows to particles to go to infinity from the disc. The 
necessary condition for this is 



GM Q.krB^ 



f 



> 0. 



(41) 



This condition means that the energy per particle is positive, 
what is necessary ( but not sufficient) to have positive v'^ 
at large distance from the source. The substitution of the 
explicit dependance of M into Eq.(Uoh gives that 



E : 



{2a^ - 3) 



GM 
2r 



(42) 



and the outflow of the wind at large distance from the central 
source is possible only under the condition a > -^3/2. 

The integral relationships between the disc and the 
wind at the base allows us to define the velocity of plasma in 
the wind at the base. Transformation of equation ( po| ) gives 
that on the disc surface 



0. 



(43) 



dr 2 ^ dr 2 ^ dr 

Here we use the fact that = 

Let us multiply equation ( ^sj ) by Q.k and taking into 
account that rQ.k = Vk transform the equation to the form 



GM 



dM ^ 

'~a — ' 

dr 



-M- 



j^d_Yk_^dM_ 
dr 2 dr 



VkV. 



VkB, 
f 



0.(44) 



8V,f 



If to take into account that — 

r 

above can be transformed as follows 

The substraction of Eq.(^) from (^) gives 
dM Vf dM ( 



^Q-, the equation 
VkB^ 



f 



= 0. 



(45) 



dr 2 



+ 



dr 



Vkv^ = 



(46) 



2 2 I 



The velocity of the wind at the disc surface v 
Therefore, after simple transformations we obtain from 
that 



Vp + (% - Vkf = 0. (47) 

It follows from (^) that at the base of the wind from the 
disc the poloidal velocity of the plasma equals to 0, and the 
azimuthal velocity equals to Keplerian velocity of rotation. 
Physically this means that the wind fully corotates with the 
disc at the base. 



4 SELFSIMILAR SOLUTIONS 

The accretion rate M appears power law function of the disc 
radius if the Alfvenic radius is proportional to the radius 
where the field line is rooted in the disc. Therefore, it is 
interesting to look for selfsimilar solutions of the problem 

We are primarily concerned here with the plasma dy- 
namics immediately above the disc, at distances z much 
smaller than the radius of the disc -Rdisc. In the limit 
z <^ -Rdisc, only two parameters with the dimensions of 
length remain in the problem: z and r. Therefore, the solu- 
tion will be self-similar in this limit (Barenblatt 1979). Let 



us underscore an important feature of these solutions: they 
describe flows only at small distances from the disc, and are 
not applicable at distances comparable to or exceeding the 
size of the disc. 



Dissipationless discs 7 



4.1 Selfsimilarity conditions 

The most comprehensive study of types of selfsimilar flows 
has been investigated by the team from Athens University 
(Vlahakis & Tsinganot 1999). Here we will be in terested in 



the selfsimilarit y of the form proposed initially by 



Blandford 



fc Payne (1982) In this kind of selfsimilarity all the variables 



depend on the coordinates in the form 



G{z,r) = 



(48) 



where z, r are the cylindrical coordinates, and S is the self- 
similarity index. 

The steady-state equations for an ideal, cold plasma 
(with pressure P = 0) can be written as 



p(vV)v = -^VBV -^(BV)B - 
on 47r 



GMH 



(49) 



This equation defines the structure of the wind. Accord- 
ing to the selfsimilarity assumption all the variables in these 
equations can be presented as 



v{r,z,(t)) = 
p{r,z) = 
B(r,z,0) = 



v{z/r,(f 
-p{z/r), 
^B{z/r, 



(50) 



This representation of the variables says that they are 
scaled as the power law of r and all functions with "depend 
on the angle ^ defined as tan^ — z/r. 

The superscripts S^, 5p, and 5b are determined from the 
following conditions. Substituting (^o|) into (^^ leads to the 
equations 



25b — 5p = 25„ = 1 . 

It follows from these equations that 

v(r, 2,0) = r~^^^ v{z/r, (j>) , 

p(r,z) = r-^p{z/r), 

_ (i+a) - 

B(r, 2,0) = 7-2 B(2/r, 0). 



(51) 



(52) 



Only one index 5 defines the family of the solutions. In our 
previous work this index has been fixed by the condition 
( p^ and (O) . Formally it follows from these conditions that 
vB ~ R~ and therefore 5 = 0. However, it is worth to 
pay attention that in the relationship between E^p, B and 
V only the component of the velocity Vr perpendicular to 
the poloidal field line comes in. But this component does 
not play any role in the dynamics of the disc or the wind 
because we consider the case Vr <C Vk- This component of 
the velocity is present neither in the wind equations or in 
the equations connecting the disc and the wind. Therefore, 
we can relax the condition of the selfsimilarity and consider 
a more general case when the component of the velocity 
in the disc does not follow to the selfsimilarity prescriptions 
(po|). The index 5 remains free parameter of the problem in 
this case. 

Let us now consider how the index 5 is connected with 
the index a in the power law dependence of the accretion 
rate. The position of the Alfven point is defined by the con- 
dition. 



Bp 
A-Kp 



(53) 



Substitution of the conditions ( po| ) in this relationship gives 
that at the Alfven point the ratio ZA/r^ is constant for 
all field lines. All the filed lines in the selfsimilar solutions 
differ by scaling of the coordinates. The shape of a field line 
is described by the equations 



z = To g{^) tan(C) , 



(54) 



where rg is the cylindrical radius of the field line at the base. 
It follows from these relationships that at the Alfven point 
the ratio r^/rg = a is constant for all the field lines. This 
means that in the selfsimilar solutions under the considera- 
tion the accretion rate varies with radius according to ( p3|) . 

Substitution of (^3|) into ( |3^ taking into account ( |52[ ) 
gives the following relationship between a and 5 



3a' 



2(q2 - !)■ 



(55) 



Q can change from | up to oo. It is interesting that at the 
same time the index of selfsimilarity varies only in the limits 
between | and |. 

It is interesting also to estimate the dependence of radial 
velocity of the matter in the disc on radius. Taking into 
account that the toroidal electric field depends on r in the 
disc as ~ and the magnetic field varies as B ~ 

l + iS 

r 2 . The frozen-in condition gives that the radial velocity 

s~i 

depends on r as ~ r . The radial velocity falls down 
with r or remains constant at 5 = i. 



4.2 Dimensionless variables 

Is convenient to consider the flow in the dimensionless vari- 
ables. In the selfsimilar solution all the geometrical variables 
can be scaled in the units of Tq , the radius at the disc where a 
selected fleld line is rooted. The magnetic and electric fields 
are dimensionlessed by the value Bo = where ip is the 
magnetic fiux through the disc with radius rg . Kepler veloc- 
ity Vk at the radius is used to dimensionless the velocities. 
These variables provides us a characteristic density 



Po 



(56) 



which is used to dimensionless the density of the plasma. 
It follows from the assumption ^ that the dimen- 
sionless density of plasma in the disc is much greater than 1. 



5 RESULTS 

The technic of solution of the problems of selfsimilar MHD 
outflows has been develop ed in many papers starting with 
Blandford & Payne (1982). The most compre hensive study 
of the selfsimilar outflows is presented 



by 



Vlakhakis & 



Tsinganos. (1998). In these papers the methodology of the 



solution is also presented in details. 

The most difficult task at the solution of the MHD prob- 
lems of outflows is the selection of the unique solution reg- 
ular at the critical points: slow magnetosonic, Alfvenic and 
fast magnetosonic critical points. In our problem the slow 
magnetosonic critical point is absent because the outflow is 
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Figure 5. The filed lines for the solution corresponding to o = 8, 
/ = 0.02. 



2500 



2000 



1500 



1000 



500 




20 40 60 80 100 120 

r 



Figure 6. The filed lines for the solution corresponding to o = 5, 
/ = 0.1. 



cold. It is necessary to find a solution regular at the Alfvenic 
and fast magnetosonic critical points. 

Typically the selfsimilar solutions are regularized at the 
Alfvenic point. It is sufficient to regularize the solu tion at 



the Alf venic point to specify the selfsimilar solution ( Bland- 
ford & Payne 1982). The group of Athens university was 
succeeded to obtain the selfsimilar solution re gular in both 
Al fvenic and fast magnetosonic critical points ( Vlakhakis et 
al. ||200(j ). 

In our work we follow to the conventional approach and 
look for the solutions regular at the Alfvenic point only. For 
astrophysical implications the solutions with large a are the 
most interesting. We present here the results for a = 8 and 
/ = 0.02. Field lines are shown in Fig. H. The field lines start 
from the disc u nder the angle 21 



B landlord 



^ Therefore the 

& Payn^ (L982) condition of outflow of matter from the disc 
is fulfilled. Interesting feature of this solution is that initially 
the field lines of the solution diverge. However, at some dis- 
tance from the disc they start to recoUimate forming the 
jet-like flow. The scale along z axis is strongly reduced in 
Fig. 1^. The solution reaches distances greater than million 
times of the initial radius. Nevertheless, the flow does not 
cross the fast mode critical point. 

This behavior is typical for large a. For smaller a the 
field lines looks like in the solutions obtained by Blandford 
& Payne (1982). The field lines for a = 5 are shown in Fig. ^ 
They start from the disc under the angle 16° . 

Such dependance on a of the solution is understandable. 
The larger a the larger the magnetic field or smaller the 
mass flux density along a fleld line. In any way this results 
into stronger affect of the magnetic field on the fiow. As 
the result, at larger a the effect of the magnetic coUimation 
appears stronger. 

In Fig. we show the variation of the plasma velocity 
along a field line with z. The growth of the velocity occurs 




Figure 7. Variation of velocity with z for a = 8, / = 0.02. 
Solid line - full velocity, dashed line - poloidal velocity, dotted 
line - azimuthal velocity. Horizontal dashed-dotted line shows the 
upper limit on the velocity following from Eg. (^) 



at z comparable with r. It is interesting that the terminal 
velocity more than 10 times exceeds the Keplerian velocity. 

Surprisingly, the transformation of the Pointing flux 
into kinetic energy flux in the appeared very efficient. The 
ratio of the Pointing flux over the kinetic energy flux in the 
winds conventionally characterized by the a - parameter. 
The behavior of this parameter along a field line is shown 
in Fig. ^. Initially the wind is strongly magnetically domi- 
nated, but at high distances the energy is fully transformed 
into bulk motion kinetic energy. 



Dissipationless discs 9 




Figure 8. Dependance of the magnetization of plasma on z. 



6 DISCUSSION 
6.1 Energy budget 

According to our solution the accretion rate M goes to zero 
at the gravitation centre. This means that no matter falls 
onto the centre. All the accreted matter is ejected. Never- 
theless, in the solutions with > 1.5, energy per particles 
in the wind is positive and therefore the wind carries pos- 
itive energy. This produces impression that conservation of 
the energy is violated. 

Apparent energy conservation violation certainly de- 
serves of a detailed consideration. First of all we have to 
point out that the violation of the energy conservation in 
our solution is impossible, because the energy conservation 
law has been used explicitly in the solution. In the wind the 
energy is conserved by imposing that E from (20|) is con- 
stant along any field line. At the same time the dynamics of 
the wind is connected with the dynamics of the disc by the 
conservation laws: mass, angular momentum and energy. 

The conservation of the energy in the disc - wind system 
follows from equation (^). Let us select some fragment of 
the disc, surrounded by the surfaces Si, S2 and S3 as it is 
shown in Fig. ^. According to ( ^o| ) the energetic budget of 
this fragment is as follows 

M{ri) ^ - M{r2) ^ + {Epv.SaU^a = 0. (57) 
2ri 2r2 

This relationship is valid for the disc and the wind at 
the base. Nevertheless, the values E and pvpS-j are con- 
served along filed lines. Here Vp is the poloidal velocity of the 
wind and S'3 is the square between two neighbour field lines. 
Therefore, in Eq.(p7|) we can replace EpvzSa on the value 
EpvpS-i which correspond to any point located arbitrary far 
from the disc. Equation ^t\) can be summed over all similar 
fragments of the disc and attached wind starting with some 
value of r2 up to the outer edge of the disc. After that only 
two terms along surfaces S2 and W survive. Integrals along 
surfaces A and B equal to zero because there is no matter 
flux through them. Integral along the surface So located at 
the outer edge the disc equals to zero because = on this 
surface. Thus we have 




Mir2) — + 



Epvp dS^O. 



(58) 



Figure 9. The integration of the energy flux over the surface sur- 
rounding the disc. The surface is shown by thick lines. The gas 
inflow into the disc at it outer edge with the energy per particle 
-5 = gives zero contribution into the total budged. E > ev- 
erywhere in the wind. However, the matter crossing the surface 
S2 carries flux of negative energy which totally compensate the 
flux of energy of the wind. 

The value Epvp dS = Lw is the flux of the energy in 
the wind which is positive in the case under consideration. 
It is exactly balanced by the flux of negative energy which 
flow to the gravitation centre through the surface 82- Thus, 
total flux through any surface containing the disc and the 
wind as it is shown in Fig. ^ equals to zero. This means that 
all the energy flux in the wind is the gravitational energy 
released at the accretion of the mass remaining in the disc. 

The solution of the paradox with the apparent violation 
of the energy conservation directly follows from Eq.(p8[). The 
mass flux onto gravitation centre M(r) goes to zero at the 
centre indeed. Nevertheless, the negative energy fiux which 
is the product —M{r)^^ diverges at the centre for > 1.5. 
This flux of negative energy exactly equals the positive flux 
of the energy of the wind. 

The budget consideration demonstrates that the solu- 
tion describes the process of transfer of the energy of matter 
infalling onto gravitating centre to the matter ejected from 
the disc. The magnetic field provides the mechanism for this 
transfer. 

This paradox physically is illustrated by Fig. |l^. In 
some sense it is analog of well known Penrose process of 
energy extraction from a black hole (Penrose, 1969). Let 
us imagine that some amount of matter is dropped onto a 
gravitating centre with zero initial energy. On the way to 
the gravitating centre infinitesimally small amount of mat- 
ter (say one hydrogen atom) splits from the matter and falls 
down onto the centre separately. The energy of this atom is 
Efj — —^YF S'lid goes to minus infinity with r 0. The 
energy equal to —E^ is transferred to the remaining mass by 
means of the magnetic field. Due to this the remaining mass 
will be ejected from the gravitation centre with the energy 
—Eh- While the hydrogen atom approach to the gravitation 
centre. Eh goes to —00. This means that the remaining mass 
will be ejected with infinite positive energy. 

This way infinitesimally small amount of mass falling 
onto the gravitating centre results into ejection of finite mass 
with infinite energy. Of course this is the result of gravita- 
tional singularity. In application to the real astrophysical ob- 
jects it is necessary to take into account that some minimal 
possible radius always exists. In the case of stars (neutron or 
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Figure 10. Illustration of energy conservation at tiie ejection- 
accretion processes. The matter shown by shadowed region falls 
initially with zero energy. On the way to the gravitation centre 
one hydrogen atom (infinitesimally small part) splits from the 
matter and falls separately. Magnetic field provides mechanism 
of the e nergv transfer from the atom to the rest of mass which 



is ejecte ^with infinitely large energy while the atom approaches" 
the centre. 



protostars) this is the radius of the stars, in the case of black 
holes this is Schwarzschild radius. If so, the accretion disc 
has some inner radius ri„ and the relationship between the 
energy flux in the wind and the accretion rate is as follows 



Lw = M 



GM 



■ V? 

= 



The relationship ( |42D gives EparUcie in the wind, 
useful to present this energy in the form 



Eparticle — (2q 3) 



(59) 



It is 



(60) 



There is rather conventional opinion (Livio) that the 
energy of the particles in the wind from the disc should be of 
the order of Keplerian velocity. It follows from (^) that the 
particle energy can remarkably exceed the Keplerian energy 
provided that well exceeds 1.5. 

6.2 Efficiency. 

The efficiency of the outflow formation by a system can be 
deflned as the ratio of the kinetic luminosity of the outflow 
( jet plus wind) over the energy released at the accretion. 
According to Eq.(|5^) the efficiency of the system under the 
consideration is equal to 1. Almost 100% of accretion energy 
is transformed into the energy of the outflow. 

At the estimate of the efficiency of real objects bolomet- 
ric luminosity is used as the estimate of the energy released 
at the accretion. This estimate follows from shakura (1972) 
solution. Therefore, the efficiency rj equals to 

(61) 

In our solution bolometric luminosity equals to zero be- 
cause all the energy goes to the energy of the wind and 
no processes which transform the energy of plasma into the 
thermal energy. Formally, for our solution -q — <x. However, 
in real objects dissipative processes, viscosity and finite elec- 
tric conductivity, always takes place and some fraction of 
energy will go into the radiation of the system. But the en- 
ergetics of this radiation does not already directly connected 



V 




Figure 11. Distribution of column density along line of sight for 
the solution with a = 8, / = 0.02. 



with the en ergetics of accretion like in the Shakura & Sun 



yaev (1973) solution. Therefore, we can expect that for some 



astrophysical objects 77 could be very high. This conclusion 
naturally explains the ratio of the kinetic luminosity on the 
jet from M87 over the bolometric luminosity of this object 
which archives the value as much as 100. 



6.3 Jet formation. 

The case when a well exceeds -^3/2 is especially interesting 
for astrophysical implications. Such a wind from the accre- 
tion disc provides the rate of accretion practically constant 
along the radius of the disc. Indeed, the disc wind with a = 5 
gives that the rate of accretion varies as M oc r^''^**, while for 
Q = 8 the dependance becomes M oc r^^^^^. The accretion 
rate appears almost constant in the largest part of the disc 
and sharply reduces at the centre. This means that the out- 
flow rate is weak from the disc and the most dense and fast 
outflow occurs from the very central part of the accretion 
disc. The matter flux density varies as 



pv = 



2ic 



1 



\ I max / 



(62) 



It follows from this equation that increase of a reduces the 
mass flux from the disc. This is natural, because increase 
of a means increase of the angular momentum carried out 
by every particle of the wind. Therefore, the same angular 
momentum can be carried out by less amount of the wind. 

Obtained form of the accretion rate explains rather puz- 
zling properties of the jets. Jets always ar formed at the cen- 
tral part of the system producing impression that some way 
the central source is responsible for the jet production. In 
the solution we see that practically all the outflow concen- 
trates at the very central part of the disc. 

To illustrate this property of the solution we have calcu- 
lated for the solution with a = 8 the distribution of column 
density along line of sight. The viewing angle of the disc is 
42° . Integration of the density near the disc was limited by 
the points on the fleld lines which are located at the Alfvenic 
radiuses. Therefore the disc appeared visible only at the cen- 
tral part. Otherwise integrals diverges at the disc. The result 
of calculation is presented in Fig. hll This figure well illus- 
trates that the obtained solution naturally explains why the 
jets are formed at the very central part of the accretion disc. 
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7 CONCLUSION. 

The selfconsistcnt solution of the problem of the disc ac- 
cretion and plasma outflow from the disc shows that such 
puzzling properties observed in real astrophysical objects 
as high energetic efficiency of jets, launching of jets from 
very central part of the accretion discs and high velocities of 
outflows which in some objects well exceeds the Keplerian 
velocities, are naturally explained in frameworks of general 
approach assuming that the largest part of the angular mo- 
mentum from the disc is carried out by the wind rather than 
due to viscous stresses. Although a significant work should 
still be done to explore the astrophysical implications of the 
obtained solution to specific astrophysical objects and gen- 
eralization of the solution for relativistic objects, we believe 
that this approach is rather promising for explanation of 
properties of jets from objects of different nature, starting 
with jets from protostars up to jets from AGN. 
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